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We show how the effective equation of motion for a distinguished quantum particle in an ideal 
gas environment can be obtained by means of the monitoring approach introduced in [EPL 77, 
50007 (2007)]. The resulting Lindblad master equation accounts for the quantum effects of the 
scattering dynamics in a non-perturbative fashion and it describes decoherence and dissipation in 
a unified framework. It incorporates various established equations as limiting cases and reduces to 
the classical linear Boltzmann equation once the state is diagonal in momentum. 
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I. INTRODUCTION 

A basic problem in the field of open quantum dynam- 
ics is the question how the motion of a tracer particle, 
such as a Brownian particle, is affected by the presence 
of a background gas [l| . More specifically, one may con- 
sider a single distinguished test particle which moves in 
the absence of external forces, but is interacting with an 
ideal, non-degenerate, and stationary gas. The elastic 
collisions with the gas particles will affect the motional 
state of the tracer particle, and we are interested in the 
appropriate effective equation of motion for its (reduced) 
density operator which incorporates the interaction pro- 
cess in a non-perturbative manner. This master equation 
is necessarily linear, since it pertains to a single particle, 
and it is aptly called, in analogy to the case of a classical 
tracer particle 0], the quantum linear Boltzmann equa- 
tion (QLBE). However, one should not confuse it with 
a linearized quantum equation for the single particle gas 
state of a self-interacting quantum gas, sometimes called 
by the same name (though the notation "linearized quan- 
tum Boltzmann equation" would seem more fitting). 

The dynamics to be described by the QLBE can be 
quite involved because the tracer particle may be in a 
very non-trivial motional state, characterized for example 
by the non-classical correlations between different posi- 
tion and momentum components found in a matter wave 
interferometer On the long run, the tracer particle 
will approach a stationary "thermalized" state, while the 
ever increasing entanglement with the gas will reduce its 
quantum coherences already on much shorter time scales. 
A limiting case occurs if the tracer particle can be taken 
as infinitely massive, so that energy exchange during the 
collisions can be safely neglected. In this case one expects 
pure coUisional decoherence, i.e., a spatial "localization" 
of an extended coherent matter wave into a mixture with 
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reduced spatial coherence. 

This problem was first investigated by Joos and Zeh in 
a linearized description However, a non-perturbative 
treatment is required to describe how the spatial coher- 
ences in an interfering state get reduced the more the 
better the scattered gas particles can "resolve" the dif- 
ferent interference paths, and to account for the satura- 
tion of this effect with increasing path difference 0] . 
This loss of coherence, which may be related to the 
"which path" information revealed to the environment, 
was observed experimentally with interfering fuUerene 
molecules in good quantitative agreement with decoher- 
ence theory [7,]. 

The situation is much more involved if the ratio m/M 
between the mass m of the gas particles and the mass 
M of the tracer particle cannot be neglected. In this 
case the particle experiences friction, it will dissipate its 
energy and finally thermalize. The appropriate effective 
equation must then be able to describe the full interplay 
of decohering and dissipative dynamics. An important 
advancement in this direction was the proposal by Diosi 
8] of an equation based on a combination of scatter- 
ing theory and heuristic arguments. In this derivation a 
number of ad-hoc approximations had to be introduced 
when incorporating the Markov assumption in order to 
end up with a time-local master equation in Lindblad 
form. As is notorious in non-perturbative derivations of 
Markovian master equations, these approximations are 
not unambiguous and very hard to motivate microscopi- 
cally. 

One way to overcome this ambiguity problem was re- 
cently proposed by one of us [9j. This method, called the 
monitoring approach, treats the Markov assumption not 
as an approximation to be performed when tracing out 
the environmental degrees of freedom, but incorporates 
it before this trace is done by combining concepts from 
the theory of generalized and continuous measurements 
with time dependent scattering theory. When applied to 
the present case, the essential premise of this approach 
is to assume that both the rate and the effect of indi- 
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vidual collisions between the tracer and the gas parti- 
cle are separately well defined. The Markov assumption 
then enters by saying that three-particle collisions are 
sufficiently unlikely to be safely neglected, as are subse- 
quent collisions with the same gas molecule within the 
relevant time scale. This assumption excludes liquefied 
or strongly self-interacting "gas" environments, but it 
seems natural in the case of an ideal gas in a station- 
ary state. The only real freedom in this framework of 
the monitoring approach lies in the choice of two micro- 
scopic operators. Selecting the operators suggested by 
microscopic scattering theory then leads to the equation 
in an unambiguous way. 

The present result was already announced in [lo| . Here 
we give a more detailed derivation^ , presenting two inde- 
pendent ways of evaluating the environmental trace. We 
will also point out that various limits reduce the QLBE 
to well-established results. In particular, one obtains the 
weak-couplin g ve rsion of the QLBE, proposed earlier by 
one of us pH Il2l. [l3j , if the appropriate limit is taken by 
replacing the scattering amplitudes with their Born ap- 
proximation. Other limits lead to the generalized form 
of the Caldeira-Leggett master equation , the master 
equation of pure coUisional decoherence, and the classical 
linear Boltzmann equation. 

The structure of the article is as follows. In Sect. Ullwe 
briefly review the monitoring approach and specify the 
microscopic operators for the problem at hand. Before 
delving into the calculations we present the form of the 
resulting QLBE in momentum representation in Scct. lIIll 
This allows us to discuss the relation of the QLBE to 
the classical linear Boltzmann equation. Section HVl then 
starts out with the calculation in momentum basis and 
explains why a straightforward evaluation of the trace is 
impossible. A first remedy, based on the restriction to 
wave packet states of incoming type is given in Sect. [V] 
Section IVII provides an alternative way of doing the en- 
vironmental trace, which is based on a formal redefini- 
tion of the scattering operator. Section IVIII is devoted 
to calculating the coherent modification part of the mas- 
ter equation using the same wave packet technique as 
in Sect. [V] The basis independent "operator form" of 
the QLBE is obtained in Sect. IVIIII it shows immedi- 
ately that the master equation provides the generator of 
a completely positive and translationally invariant quan- 
tum dynamical semigroup. Section IIXI summarizes the 
various limiting forms of the QLBE, and we present our 
conclusions in Sect. 1x1 



^ We emphasize that the word "derivation" is used here in the 
physicist's sense, implying that arguments and approximations 
are invoked which-though physically stringent and leading to a 
uniquely distinguished equation-may be hard to substantiate in 
a proper mathematical framework. We certainly do not claim 
to provide a mathematically rigorous proof, noting that even 
the classical Boltzmann equation still lacks such a mathematical 
derivation. 



II. THE MONITORING MASTER EQUATION 
FOR A TRACER PARTICLE IN A GAS 



A. The monitoring master equation 

Let us start with a brief review of the monitoring ap- 
proach ^9(1 . It yields a Markovian master equation that is 
specified, apart form the system Hamiltonian H, in terms 
of two operators, a rate operator f and a scattering op- 
erator S. 

The operator f is positive and in the present context 
it has the defining property that its expectation value 
yields the probability of collision with the gas particles 
in the small time interval At, 

Pr(CAt|P'»Pga.) = Tr(r[p®pgas])Ai + 0(At2)(l) 

Here, p is the system density operator which describes, in 
the present application, the motional state of the tracer 
particle. The operator pgas is the effective single parti- 
cle state of the gas environment, and it is assumed to be 
stationary (but not necessarily in thermal equilibrium). 
Thus, r acts in a two-particle Hilbert space, and its task 
is to incorporate the tracer state-dependence of the col- 
lision probability into the dynamic formulation. 

The scattering operator S, on the other hand, is uni- 
tary, and by definition it yields the two particle state 
after a single collision, so that, upon tracing over the gas 
particle, we obtain the new tracer particle state (in inter- 
action picture) after a single scattering event took place 

[a, 

p' = Trgas(S[p®Pgas]St). (2) 

The monitoring approach Q now implements the Markov 
assumption by combining the state dependence of the 
collision probability ([T]) with the transformation ^ in 
a way which is consistent with the state transformation 
rules of quantum mechanics [l6| . using concepts of the 
theory of generalized and continuous measurements [13, 
[isl . figj l . In the Schrodinger picture one thus obtains the 
effective equation of motion 

^p = ^[H,p]+Cp + np. (3) 
at in 

The superoperators C and TZ are best specified in terms of 
the nontrivial part T of the scattering operator S — l-|-zT. 
The part Cp then takes the form [S] 

Cp = Trgas [p <g> pga,] fi/^jt) 

-^Trga3(r/2TtTrV2 [p^pgj) 
-^Trga3([p®Pgas]ri/'TtTr/2y (4) 

It describes the incoherent evolution of p due to the pres- 
ence of the gas environment. The part TZp, on the other 
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hand, is given by^ 



7^p - iTrgas([ri/'Re(T)ri/ 



gas 



(5) 



where Re (T) = (T + T^) /2. It is responsible for a uni- 
tary modification of the evolution, a renormalization of 
the system energy due to the coupling with the environ- 
ment. 

We would like to emphasize that the evolution de- 
scribed by ([3]) is non-perturbative in the sense that the 
collisional transformation described by S is not assumed 
to be weak. Moreover, note that the incoherent part (g]) 
is manifestly Markovian even before the environmental 
trace is done. 



B. Rate and scattering operators 

In the framework of the monitoring approach the only 
essential freedom lies in the choice of the operators f, 
T, and H appearing in Eqs. ([S])-©. In this section we 
will specify them on a microscopic basis. Before that it 
is helpful to consider with some care Pgas, the effective 
single particle state of an ideal gas with number density 

To be describable by a normalizable state, the gas must 
be confined, say with periodic boundary conditions, to a 
finite spatial region O with (large) normalization volume 
Let us denote the projector to this spatial region as 



Ic 



da; la;) lx\ 



(6) 



Using the double- bracket notation for the volume- 
normalized momentum states, the density operators cor- 
responding to these proper vectors take the form 



Pp 



\p))M 



\n\p){p\\ 



(7) 



X a n. The most natural choice for /i is of course the 
Maxwell-Boltzmann distribution, see ([M]) below, but we 
will keep /j, unspecified in order to indicate that the par- 
ticular form of /i is not relevant for most of what follows. 

In principle, projections similar to the Iji in ([8]) are also 
needed when defining the operators f, T, and H of ©-(O. 
To avoid clumsy notation we will instead present them in 
their unrestricted form and take care of the restrictions 
during the calculations below. 

Since the tracer particle is supposed to move in the 
absence of external forces the Hamiltonian part of ^ is 
given by H = P^/2M, where P is the momentum operator 
of the tracer particle. The two-particle operators V and 
T depend on the relative coordinates between tracer and 
gas particle, and it will be convenient to denote relative 
momenta by 



rel(p,P) := 



-P 



M' 



(9) 



with = TfiMj (M + m) the reduced mass. Thus, the 
momentum dyadics corresponding to the different fac- 
torizations of the total Hilbert space Htot ='H® Ti-gas = 
'Hem ® Ti-rci are related by 



\P){P'\®\p){p'[ 



\p + p){p' + p'U (10) 

®|rel(p,P))(rel {p'.P') 



In classical mechanics, the collision rate is ob- 
tained by multiplying the current density jo [Pi P) — 
'^gas Ire^Pji-*)! /m* of the relative motion with the to- 
tal scattering cross section dtot (which also depends on 
the relative momentum). It seems therefore natural to 
define f as the corresponding operator on Hem ® ^rci, 



bo (P, P) atot (rel (p, P))] 



rel 



OJrol 



with 



|rel(p,P)|atot (rel(p,P)), 



(11) 



(12) 



Here the \p) are the usual improper momentum eigen- 
vectors, {x\p) ~ {2TTh)^^^^ exp {ix ■ p/h). Since Pgas is 
stationary it must be a convex combination of the pure 
momentum states ([7]). It is completely characterized by 
the gas momentum distribution fi (p) , a positive function 
satisfying/ dp/i (p) = 1 . Thus pgas has the form 



Pgas 



J dpM(p) Pp = 



lnMp)lo, (8) 



where p is the unrestricted momentum operator of a sin- 
gle gas particle. This state is normalized, Tr (pgas) — 1, 
and it is uniform in position, (a;|pgas|a;) = for 



^ A marginally different expression was given in Ref. Q, see the 
discussion in Sect. IVIII 



Indeed, for normalized and separable particle-gas states 
the expectation value of this operator yields the collision 
rate experienced by the tracer particle, provided their rel- 
ative state is of incoming type. If the two-particle state is 
of outgoing type, on the other hand, the motion of the rel- 
ative coordinate is directed away from the origin, so that 
the particle and the gas molecule never interact. Still, 
the operator (jlip would yield a finite expectation value 
in that case, since it depends only on the modulus of the 
relative velocity and not its orientation. A proper defi- 
nition of r should therefore also include a projection to 
the subspace of truly incoming relative motional states. 
Unfortunately, it is rather difficult to formulate this pro- 
jection in a way so that one can work with it in concrete 
calculations. Therefore, instead of using a more refined 
definition we shall stick with Eq. pip keeping in mind 
that it is valid only for incoming states of the relative 
motion. 
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As the last step, we have to define the operator S = I + 
iT describing the effect of a single collision. It is natural 
to use scattering theory for a microscopic definition [3] ■ 
The center-of-mass coordinate then remains unaffected, 



[So] 



(13) 



and the scattering operator of the relative coordinates 
So = I + iTo is fully specified in terms of the complex 
scattering amplitudes / (p/ , Pili, which are determined 
by the inter-particle potential |20l |. 



f{pf,p,). (14) 



Note that a delta-function in ensures that the en- 
ergy is conserved during an elastic collision changing the 
relative momentum from to . 

The scattering amplitude also defines the cross section 
required in (fTT|) . The differential cross section is given 
by 



(^{Pf^Pt) = \f{Pf,Pt 
and the total cross section reads 



CTtot (Pj) = y 1/ P» 



(15) 



(16) 



where n is a unit vector with dn the associated solid 
angle element. 

It is important to keep in mind that the S-matrix (jl3|) 
provided by scattering theory is physically meaningful 
only for proper incoming states of the relative motion, 
even though it is defined on the whole Hilbert space. This 
is schematically shown in Fig. [T] where we contrast the 
action of So on an incoming wave packet with its effect on 
an outgoing state, showing that an outgoing wave packet 
may get spuriously transformed. The reason is that the 
M0ller operator f2_|- = limt^oo U (t) Uq {—t) used to con- 
struct So = involves a free backward evolution 
in time Uq {—t), followed by a forward motion U (t) in 
the presence of the interaction potential. To avoid this 
undesired transformation in ^ we must either ensure 
that outgoing wave packets contribute with a zero colli- 
sion rate or we have to modify So such that it leaves the 
outgoing contributions invariant. 



III. THE QLBE IN MOMENTUM 
REPRESENTATION 

Before we proceed to derive the quantum linear Boltz- 
mann equation let us present the result in the basis of 
improper momentum eigenstates \P), where it takes a 
particularly simple form. This permits to introduce the 
complex rate function Mi„ to be evaluated in Sects. [TvT - 
IVIl and to discuss its relation to the classical rate densi- 
ties of the collision kernel. 



We will see that the momentum representation of the 
incoherent part ^ of the QLBE can be written in terms 
of a single complex function, and that it takes the form 

{P\Cp\P') 

= j dQ{P~ Q\p\P' - Q)M,^ (P, P'; Q) (17) 
^{P\P^ J ^Q^j,^ {P + Q^P + Q-Q) 

_ (P^ j ^QM,^^^p^^Q^p^^Q.Qy 



The function Afin (P, P'; Q) is defined below in Eq. 
[see also (P7)) and (|39l) ]. In order to highlight the relation 
to the classical linear Boltzmann equation, let us first 
note that the master equation takes the shorter form 

{P\Cp\P') = j dQM,^{P,P'-Q){P-Q\p\P' -Q) 



[AC {P) + aC (P')] {P\p\P') 



(18) 



once we introduce 



A'C^iP) := J dQM,,,{P + Q,P + Q;Q). (19) 

As indicated by its name, this positive function gives the 
rate known from the classical linear Boltzmann equation 
for a particle with momentum P to be scattered to a 
different momentum. It involves an integration over all 
initial gas momenta Pq and all momentum exchanges Q 
subject to the restriction implied by energy conservation. 



Kit (P) 



''gas 



x6 



dpodQuipo) 



|rel(p„,P)|2-|rel(po,P) + Qr 



xa(rcl(po,P)+Q,rel(po,-P))- (20) 

Here, p is the gas momentum distribution function from 
([5]) , the function rel (p, P) is defined in , and a is 
the differential cross section dTSl) . Carrying out the Q- 
integration one can write the rate in terms of the total 
cross section 



Kit (P) 



dPoM(Po) |rel iPo,P)\ 



m, 

XCTtot (rel (Po,^')) 



(21) 



It follows that the dynamics described by the "loss term" 
in (fT7|) . (fT8|) is fully specified by the rate in the corre- 
sponding classical equation (which involves of course a 
quantum mechanical cross section). The term leads to 
a reduction of the momentum coherences {P\p\P' ^ P) 
with a rate given by the arithmetic mean of the loss rates 
of the corresponding diagonal elements, (P|p|P) and 
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FIG. 1: Action of the S-matrix when apphed to localized wave packets of the incoming and the outgoing type, (a) An incoming 
wave packet |'(/'in) is transformed in such a way that the free motion of the resulting state So|t/'in) (indicated by the dashed 
curves) converges with the dynamically scattered wave packet at large times f2(?]. (b) An outgoing wave packet, whose forward 
time evolution will be unaffected by the scattering potential, gets strongly transformed by So. This is due to the inverse time 
evolution involved in the definition of the the S-matrix. To prevent this unwanted transformation one should (i) either attribute 
a vanishing collision rate to all outgoing states or (ii) modify the So operator such that it leaves all outgoing wave packets 
unaffected. Evaluations of £p based on these two strategies are given in Sect. |V] and Sect. IVll respectively. 



{P'\p\P') , and these momentum populations, in turn, 
are depleted in the same way as the momentum distribu- 
tion function of the classical linear Boltzmann equation. 

The "gain term" in ([TT]), ([HI) is also related to the 
classical linear Boltzmann equation, but only on the di- 
agonal P = P' . As one expects, Min is positive on the 
diagonal, and equal to the rate density M^' from the clas- 
sical linear Boltzmann equation for the tracer particle to 
end up in momentum Pf after a momentum gain of Q, 

Mf,'(P/;Q) 

= M;^{Pf,Pf-Q) EE M'^'iPf-Q^Pf) 

= / dpoAi(Po) 

"1* J 

^, ^ |rel (Po - Q, Pf)f - |rcl (pp, Pf - Q)f ^ 
X a (rel (po - Q, Pf) , rcl (po, Pf - Q)) ■ (22) 



The second equality in (j22p introduces the notation 



M^' (P,; ^ Pf) Mf„' (Pf, Pf - P,) , (23) 



for the rate density corresponding to a change of mo- 
mentum Pi to Pf. It will be useful for the discussion 
of the classical linear Boltzmann equation in Sect. IIXI 
and it yields the classical out rate (|19p as M^l^^ (P) = 
J dPfM'^^{P Pf) thus ensuring the conservation of 
probability. 

For P ^ P' the function Min is in general complex- 
valued, and it has a rather complicated form when stated 
with its explicit dependence on P, P', and Q: 



Mi^{P,P':Q) = 



X/ rel Po-Q,P- 



Pii - Pn 



,rcl Po,P - 



^ P\\ P'w ] 
M 2 j 



xf* (rel Ip^-Q^P' + 



,rel Po,P' + 



Pll - P'n 



Q 



x5 



rel (^Po - Q, - rcl (^p^. 



p+p' 



J 



(24) 
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Here we used the abbreviations 

^11 - (25) 

and 

^'l - (26) 

for the contributions in P and P' paraUel to the mo- 
mentum exchange Q ^ Q. This form of M\a {P,P'',Q) 
clearly reduces to the diagonal expression (j22p when P 
approaches P' . The curious appearance of the P|| and 
P|l contributions in (|24|) ensures, in combination with 
the delta-function, that the modulus of the initial and 
the final relative momentum are equal in both scattering 
amplitudes. This will be more obvious below in Sect. II Vl 
where suitable relative coordinates are introduced. The 
energy conservation is thus manifestly guaranteed for 
each of the scattering amplitudes separately, while the 
arguments will differ in general. 

One of the important properties of the "complex 
rate" ([M]) is that it admits a factorization of the 
P- and P -dependence, which will be crucial later 
on, when we formulate the master equation in its 
representation-independent "operator form" . Specifi- 
cally, it will be shown in Sect. IVIIII that Mjn can be 
written as a two-dimensional integration over the set 
Q"*" = {p G R"^ : p ■ Q — of momenta perpendicular 
to the momentum exchange Q. This way the integrand 
factorizes into a product of P- and P'-dependent terms, 

Mi^{P,P';Q) = f dpLip,P Q;Q) 

xL*{p,P'-Q:Q). (27) 

The functions 

L{p,P-Q) (28) 
X / Uc1(p^q,P_lq) - y,rel(p^Q,P_LQ) + yj , 



involve P|| defined in p5)) and P±q := P — P\\q- 
In the representation-independent form of the master 
equation they turn into operator-valued expressions, see 
Sect. lVml 

Concerning the coherent modification of the QLBE, it 
will be shown in Sect. IVlTl that the momentum represen- 
tation of the corresponding term ^ reads 



{P\np\P') = (P) i^n (P') , 



with 



Tl I 

E^{P) = -2nh'^ / dpofiiPo) 

xRe[/(rel(po,-P),rel(po,i'))]. (30) 
This shows how the presence of the gas changes the en- 
ergy of the particle with respect to the vacuum. This 
energy shift depends on the particle momentum and is 
determined by the real part of the average forward scat- 
tering amplitude. This phenomenon is well known in the 
field of neutron and atom interference, and it is usually 
accounted for by introducing an index of refraction, see 
Sect. HXEl 



IV. EVALUATION IN THE MOMENTUM BASIS 



A. Transformation to relative coordinates 



Our main task in deriving the quantum linear Boltz- 
mann equation is to evaluate the expressions (jH) and 
([5| , which is best done in the momentum representation. 
Starting with the incoherent part Lp, the cyclicity under 
the trace yields 



{P\Cp\P') = J dQdQ' (P - Q\p\P' - Q')M (P, P'; Q, Q') 

-\ J dPo {Po\p\P') J dPf M{Pf,Pf-Pf^P^,Pf-P) 
-\ J dP', {P\p\P'o) J dPf M {Pf, Pf - Pf - P', Pf - P',) 

with 

M (P, P'; g, Q') = (P| Trg,, {jV^'^ [\P - Q) (P' - Q'\ <E> pg,,] rV^jt) |P'). (31) 
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Upon inserting the stationary gas state ^ into (pij) we 
can simplify the expression by transforming from the two- 
particle coordinates to the center-of-mass and relative 
coordinates using pH)) . Since f and T depend only on 
the relative motion, sec pT|) and p^ . one thus finds 



M{P,P';Q,Q') 



= 6{Q-Q 



{2TThy 



\n\ 



dPo M (Po) 



X (rcl (po - Q, P) ITory Vel (Po, J' - Q)) 
X (rel (po, P' - Q) ir'J'Tl] rel (po - Q, P')> 



J(Q-Q')Mi„(P,P';Q), 



as anticipated above in (|T7| . 

It is now helpful to introduce functions of Pq , 



Pf = i-ei(po-g, — ^ — 



(32) 

(33) 
(34) 



which denote the mean of the pairs of initial and final 
relative momenta appearing in ([5^ . We also set 



q = rcl 0, 



P P' 



(35) 

These definitions imply the relations 

Pf + q = rel(po - Q,P) 
Pf-q = rcl(pQ - Q,P') 
Pi + q = rel(po,P-Q) 
Pi-q = rcl(po,P'-Q) 
Pi~Pf = Q> (36) 

which are noted here for later reference. Moreover, for 
given q we shall write 



Q-[Pf- P 



- [Pf - P^ 



(37) 



{pf~p^y 

Q± = Q-Q\\ (38) 

to denote the components parallel and perpendicular to 
the momentum exchange Q = Pi — Pj- 

The complex rate density defined in the second equal- 
ity of ((32l) now takes the form 

i2Trhf 



M,^{P,P'-Q) = J dpo^(po 



\n\ 

1/2, 



:{Pf + q\TQr„' \p, + q 



x(p,;-g|ry^Tj|p^-q'). 

We can write it as the average over the gas momentum 
distribution function of a rate density in the center-of- 
mass frame, 

Mi„(P,P';Q) = ydpo/i(Po)min(ppP,;g),(39) 



thus formally introducing 



[Pf,P^■,Q) 



(P/ +g|Tory^|p, + q> 
x(p.-g|rJ/'Tt|p^.-q> (40) 

1/2 

in terms of matrix elements of ToFq with respect to the 
relative momentum coordinates. This expression should 
be viewed here as a generalized function in the sense of 
distributions, with independent variables Pf,Pi, and q. 

The main aim of the following sections is to show, in 
two independent lines of argument, that the expression 
(l40l) should be understood as 



mn,{Pf,P^,q) -<| Q otherwise, 

(41) 



with 

mn {Pf,Pi;qA 



''gas 



m 

X 



I I / [Pf + 1^,Pt + 9±J 



/* (p/-q±,p»-qj 



(42) 



Note that this term involves a single delta-function and 
the abbreviation q^ just defined in ([38|l . 



B. Diagonal evaluation of the trace 



As a first step, let us try to evaluate (|40)) in a straight- 
forward fashion by maintaining that the operator Fq is 
diagonal in the relative momentum coordinates. Equa- 
tion then implies 



vl%) = v/Mp)Ip) 



with 



To (P) 



IpI c^tot (p) • 



(43) 



(44) 



Noting that the Tq matrix elements are given by one 
thus obtains the generalized function 

min {pj,Pi]q)^ (p^ + q) rj/^ (p., - q) 

X /{Pf+Q, Pi + q) f* {Pf - Q, Pt - q) 
2nh (p]-pI , . \ 




Clearly, the two delta-functions ensure that the en- 
ergy is conserved in each of the "elastic collision tra- 
jectories" expressed by the arguments of the two scat- 
tering amplitudes in (|45p . Employing the relation 
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5 {a + b/2) 5 {a — b/2) = S {a) 5 (6) we obtain the equiva- 
lent form 

Win {Pf,P^; q) ( ^^ ^^" ^ r^^ + q) r^^ (p^ - q) 

^ f{Pf + Q,Pi + q) f* {Pf Q,Pi q) 
xS{{pf-p.,)-q)^^. (46) 

The first delta function now requires Pj and p^ to have 
equal length. These are the mean relative momenta of 
the pairs of scattering trajectories, as defined in Eqs. ([55]) 
and ([34)1 . Given |pj — \pf\, the second delta function 
ensures that the energy is conserved in each of the scat- 
tering amplitudes individually, by granting that q, which 
expresses a "distance" between the two pairs of scatter- 
ing trajectories, is orthogonal to the momentum exchange 
Pj — Pf- The fact that possible parallel components of 
q — q± + cannot contribute to an integral over the 
generalized function (|46p can be made manifest by re- 
placing the q's outside of the delta function by the or- 
thogonal component defined in ([55]) . In other words, 
the statement P5|) is tantamount to 

Win iPf,Pt\q) 

= 5(^^^yrip,+q,)Trip,-q,) 

X f {Pf+ Q±,Pr + Qa_) f* {Pf - Q±,Pr - Q±) 

x-5((p,-p,).g)^, (47) 

which implies that the expression vanishes for <7|| 7^ 0, 
as stated in (HH). In fact, even if the integration over 
TOin involves a smooth function g (q) the second delta 
function will enforce that the latter contributes only with 
the orthogonal component of q, 



{pf,pc,q) g{q) 



[Pf,Pt;qA 



'9(g±)-(48) 



One observes on the right hand side of (|T7)) that already 
the first two lines now manifestly ensure the energy con- 
servation of the pair of collision trajectories described by 
the two scattering amplitudes. This is a crucial require- 
ment since the scattering amplitudes are not defined off 
the energy shell (notwithstanding the fact that analytic 
continuations are often considered and helpful in scat- 
tering theory). At the same time this implies that the 
physical relevance of the second delta-function has been 
accounted for once the parallel components of q have 
been set to zero. Hence, the third line in (|47p is essen- 
tially dispensable, which is all the more important since 
it renders min an ill-defined expression due to the ap- 
pearance of the arbitrarily large normalization volume 

As is well understood, the evaluation carried out in this 
subsection does not yield a well-behaved result because it 



takes the momentum-diagonal form (|12p of the rate op- 
erator r too seriously. It was already discussed in Sect. 
Ill B I that either f should involve a projection to the sub- 
space of incoming wave packets, or that the operator S 
should be redefined such that it keeps the outgoing wave 
packets unchanged. These two strategies will be imple- 
mented in Sects. IVl and I VII yielding identical results. As 
one expects, the overall structure of (|17|) . which is dic- 
tated by the energy conservation, will not change, but 
the third line will be replaced by a proper normalization. 



V. WAVE PACKET EVALUATION 

The aim of this section is to evaluate the generalized 
function 



{pf,p^■,q±) 



{2TTh)- 



Ap^-qA^y^A\Pf^qA} m 



by consistently incorporating the fact that the rate op- 
erator To should have a vanishing expectation value for 
those states of the relative motion that are not of in- 
coming type. As a first step, we will write as the 
expectation value of a non-hermitian operator with re- 
spect to a properly normalized momentum state of the 
relative motion. For that purpose it is convenient to in- 

1 /2 . 

troduce the operator Zq ToFg and its translation by 
the momentum q^. 



Zq = cxp ( -i 



oFq' exp [i — J (50) 



T, 



where x^oi is the position operator of the relative coordi- 
nate. Moreover, we note that, analogous to ([71), a volume- 
normalized momentum state of the relative motion has 
the form 



Pp. 



\n\Pi){p^\\s 



(51) 



Combining ([50| and (|5ip one finds that the complex rate 
density (|49p can be taken as the diagonal momentum ma- 
trix element of an operator product, min {Pf,Piiq±) = 

(pj|Zq^/9p.zLq^ \Pf), provided the projection to the nor- 
malization volume is included. If we further denote the 
projector to improper momentum eigenstates as Pp^ = 
\pj){Pf\ we can write 



mi 



{Pf,P^■, 



Trfzt 



(52) 



The complex rate density has now the form of an 
expectation value with respect to a state pp_ which is 
properly normalized, Tr (pp.) = 1. As discussed in Sect. 
IIIBl the rate operator Fq, and therefore also the non- 
hermitian operator Z^.^^ Pp^ Zg^ should include a restric- 
tion to the subspace of truly incoming relative motional 
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FIG. 2: The incoming relative momentum p^ defines a cylin- 
der with base area Ep. and height Ap.. This spatial region is 
used to implement the phase space restriction of the Wigner 
function to incoming states. 



states. Starting from (|52p . one can now implement this 
restriction in a rather transparent and intuitive fashion 
by considering the phase space representation of pp., as 
shown in the next subsection. A similar method was al- 
ready successfully applied in Q, where the effect of a 
gas on the internal dynamics of an immobile system was 
discussed by combining the monitoring approach with 
scattering theory. 



A. Phase space restriction to incoming wave 
packets 

The operator pp in ([5^ characterizes the motional 
state of the relative coordinates between particle and gas 
prior to a collision. According to (|5ip it is given by a 
plane wave which extends through the whole normaliza- 
tion volume, and it is therefore clearly not of the incom- 
ing type required for the application of scattering theory. 

The Wig ner- We y l formulation of quantum mechan- 
ics [U, |22|. [23I [24 | suggests a way to treat this prob- 
lem. Continuous variable states may be represented by 
the phase space quasi-probability function Wp {x,p) :— 
{2Trhy^ Jdqex.p{-iq ■ x/h) {p-q/2\p\p+q/2). For the 
state (|51l) the associated Wigner function reads'^ 



Wp^ {x,p) 



Xn jx) 



Sip- Pi 



(53) 



where xn is the characteristic function of the normaliza- 
tion volume \n\. Given that expectation values like ([5^ 
may now be calculated as phase space integrals, it is nat- 
ural to implement the restriction by confining (|53p to the 



^ This follows with the approximation XQ ~ f ) Xn (^ + f ) — 
xn xn (^)i which is permissible since the normalization region 
Q may be taken arbitrarily large. 



phase space area of incoming wave packets, i.e.. 



IS, 



Here, the product XA^. {^\\Pi) XSp (^-Lpi)is the charac- 
teristic function of a cylinder pointing towards the ori- 
gin, see Fig. [21 It describes those points x = a;||p. -I- x±p. 
in position space which will pass the vicinity of the ori- 
gin when propagated in the direction given by p^. Ep. 
is the base surface of the cylinder and its area will be 
taken to be equal to the total cross section below, i.e., 
|I]p. I ~ o'(pj). The interval Ap. specifies the cylinder 
height; its precise p.j-dependence will drop out of the cal- 
culation later on. 

The operator corresponding to (I54p reads 



da; 



Pi 



dx 



dw 



'Ap. I^Pil J^v, I^P,I 

^exp{i^—j\Pi-—){p, + —\, (55) 



so that compared to the unrestricted expression corre- 
sponding to ([55]) 



da; 

nW\ 



dw exp 



fx ■ w\ 



,(56) 



the spatial average over the whole normalization volume 
is simply replaced by an average over the cylinder, and 
the norm is indeed preserved, 



da; 



1. (57) 



It should be noted, though, that strictly speaking nei- 
ther Pp. nor Pp are legitimate quantum states since they 
combine a precise momentum with a finite position vari- 
ance. They should rather be seen as convenient basis 
states admitting to average over the momentum distri- 
bution function, see Eqs. ^ and ([55)1 . 



B. Restricted evaluation 

Inserting the restricted state ([55)) into ([5^ one can now 
evaluate the complex rate density to obtain a well-defined 
expression. Starting with 



[Pf,Pt;QA 



= Tr 



P 7 n' 
-q± ^Pf^<i±f^Pi 

dx 



\P^ 



dx 



-Lp. 



'Ap. I^pj JSp. \^Pi 

X {Pf + Q±\'^o^l^^\P, + Qi_ 



dw exp 



/ X ■ w\ 



-) 



w 



X (p, -q± + i^l^o^^llpf - 



(58) 



10 



we can now use with confidence the expressions (1141) and 

-1/2 



(|43p for the momentum matrix elements of To and Fq 
Thus, rTiin {Pf,Pi',Q±) takes the form 



dx 



Pi 



dec 



Pi 



dw cxp ( —i 



x ■ w 



{2TrhY 



f[Pf + Q±,P, + Q± + i:r 



X / [Pf- Q±,Pi -Q±-^ 




In the arguments of the delta function we took into 
account that (7j_ is orthogonal to the momentum ex- 
change, (pf—Pi) ■ q± = 0, as follows from Us- 
ing again the relation S {a + 6/2) S {a — 6/2) ~ S (a) S (b) 
a delta function is obtained with argument ■ w. Writ- 
ing w = w\\p^ + wj_p^, with loiip^ = {w-Pi)Pi/pf, 
this delta function renders if ||p. = 0, and as a result 
the integrand now no longer depends on a^iip.. It fol- 



lows that the integration along the cylinder axis can be 



done, J J 



We obtain 




(60) 



where we identified the cylinder base area with the to- 
tal cross section, |Sp. | = dtot (pj- One observes that the 
a; ^p. -integration over the surface Sp. of the cylinder base 
yields an approximate two-dimensional delta function in 



w±p . Combined with the one-dimensional delta func- 
tion in w\\p, = • w/pi this gives a three-dimensional 
6 (w), which permits to carry out the to-integration. We 
arrive at the well-defined expression 




It shows that the complex rate density ([5^ is essentially 
given by the product of two scattering amplitudes whose 
arguments differ in general. They correspond to scatter- 
ing "trajectories" determined by the relative momenta 
Pj, Pf, and The Pj and pj provide the arithmetic 
means of the initial and the final momenta, while Qj^ 
characterizes the distance of the "trajectories" in momen- 
tum space. Since is orthogonal to — Pj , a single 
delta-function suffices in (|61[) to ensure the conservation 
of energy in both scattering amplitudes. 

Reassuringly, this result reduces to the classical rate 
density for qj_ = 0, as can be seen easily by insert- 
ing TOjn {pf,Pi; O) into ([39)1 . This shows that the first 
line in (|6lj) may be viewed as a natural quantum gen- 
eralization of the classical case, where the "off-diagonal" 
contributions with q^ ^ represent quantum correc- 
tions. From the point of view of quantum physics, there 
is indeed no reason why the effect of the gas collisions on 
the tracer particle should be confined to the "diagonal" 
contributions given by q^^ = 0. In the present relative 
coordinate representation, the first line in (j6ip has in 
fact a straightforward interpretation. It simply provides 
the contribution of the scattering amplitudes of any pair 
of scattering trajectories, which is allowed by both the 
energy conservation and the choice of P,P', and Q in 



At the same time, one expects that a q_L-integration 
will average out the "far off-diagonal" contributions with 
large modulus |q^|, where the phases of the two scatter- 
ing amplitudes are no longer synchronous. It is there- 
fore reasonable to disregard the weak q_|_-dependence in 
the second line of ([6T|) . and this is corroborated by the 
fact that its linear dependence on q^ vanishes identically. 
This removes the second line altogether, so that we end 
up with the form claimed in Eq. (j42p . 

A noteworthy step in the present line of reasoning was 
that the base area |Sp.| of the cylinder required for dis- 
tinguishing the incoming states was identified with the 
scattering cross section a tot- This is very natural from 
a physical point of view, but it seems hard to justify 
on a formal basis. It is therefore worthwhile to present 
a second argumentation which, though very different in 
nature, leads to the identical result. 
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VI. EVALUATION WITH MODIFIED 
SCATTERING OPERATOR 

A second, more heuristic approach of evaluating TOjn 
sidesteps the issue of how to incorporate a restriction of 
the rate operator fg to the incoming wave packets, and 
takes its momentum-diagonal form (fT^ at face value. 
The unrestricted Fq then attributes a finite scattering 
rate also to outgoing wave packets, which would never 
touch the interaction region in a dynamic description. 
This forces us to consider a redefinition of the scattering 
operator, which is necessary since outgoing wave packets 
are not left invariant by the proper S-matrix Sq, as dis- 
cussed above in Fig. [TJ Let us therefore formally replace 
So by a modified operator Sq, which by construction acts 
like So when applied to asymptotic-in states, but leaves 
states with outgoing characteristics invariant. It will not 
be necessary to specify the details of this modification 
since all that is needed for the evaluation of min can be 
obtained from a single property that must hold for any 
such modified operator. It is the isometry of Sq with re- 
spect to the set of volume-normalized momentum states. 

The advantage of replacing So by Sq is that we can now 
use plane waves not merely as a basis, but as represent- 
ing proper states, because the unwanted transformation 
of their outgoing components is now formally excluded. 
Like in Sect. [TTl we use double brackets to denote mo- 
mentum states which are normalized with respect to the 
volume f2 and subject to periodic boundary conditions 
on its border. Due to the finite size of fl they form a 
discrete basis : p G P}, decomposing the identity 



as 



(62) 



Heuristically, one may view each p g P as labeling a dis- 
tinct lead connected into and out of the scattering center. 
An important difference with respect to a continuous de- 
scription is that the unitarity of the proper S-matrix, 
which expresses itself in the optical theorem, cannot be 
accommodated within this discrete setting. The optical 
theorem quantifies the diffraction limitation of the scat- 
tering cross section, telling 'how much' of a plane wave 
would pass the scattering center without distortion. If 
we describe the scattering process in terms of the am- 
plitudes corresponding to discrete momentum states, or 
distinct leads, then the possibility of 'passing the target' 
is no longer available since any matrix element may have 
a finite amplitude. [The possibility of 'forward scatter- 
ing' \\Pi)) ~^ \ \Pi)) differs from this diffractive "passing" 
and leads to an additional phase shift, see the following 
section.] This suggests to disregard the identity operator 
in Sg, which relates to the unscattered part of the state, 
and to require of the remaining transition operator that 
it conserves the norm. 



Inserting the identity (|62p we see that the sum of the 
probabilities of scattering into the different leads equals 
1, 



E l((P|To|P. 



= 1. 



(64) 



pSPn 



This is the standard property of the transition matrix 
used to describe discrete scattering problems between a 
finite number of incoming and outgoing leads, e.g. in 
mesoscopic physics ;25;j or the field of quantum graphs 
[2^ . It seems natural to demand this relation of any 
reasonably modified operatorTp. 

The use of (l64l) is that it tells us how to normalize the 
square of Tq matrix elements with respect to improper 
momentum kets. Inspecting the momentum matrix el- 
ement of the To operator given in one finds that 
1 To P involves the square of a delta-function. The 
expression should be well-defined when using the mod- 
ified operator Tq, and the obvious choice is to assume 
that it is given by the corresponding expression with a 
single delta- function and a normalization J\f (p^) yet to 
be specified. 



i(p/inip.>p = ^ip.)s 



Pt 



\f {Pf,pt 



(65) 



Approximating the summation in (|64p by the correspond- 
ing integral one finds 

1 = / dP,^i(P,inip.)p 
= AA(pj|dp,5(^^i^^ \npf.p^)? 

= A/'(Pi) |p.>tot (Pj) ■ 

This fixes the normahzation, Af (p^) ~ [|pj|tTtot (Pi)] 
and we obtain a well-defined expression for the squared 
matrix element of the modified operator Tq, 



{2TThy 



\{Pf\T'o\p^ 



= s 



p1\ \f{Pf,Pt)\ 



C^tot {Pi) \p. 



(66) 



Arriving at this equation required a certain amount of 
heuristic argumentation. It should be emphasized that 
this expression was already used in f^, where it was 
shown to yield a localization rate for coUisional decoher- 
ence that is equal to a wave packet calculation similar to 
the one in Sect. |Vl 

If we accept ((66)) the evaluation of the complex rate 
density rriin can be done in a rather straightforward fash- 
ion. Using the unrestricted rate operator To and the mod- 
ified Tq instead of To , the complex rate density from ([40]) 
takes the form 



Pf,Px,Q) ^ ^l^^{P^+Q)'^0^{P^~Q)^{Pf,P^lQ 



\S'o\\p^))\ 



\T'o\\P^))\\ 



1/2 



1. 



(63) 



(67) 
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with the formal expression 



-(P/ + g|To|Pi + q) 
x(P/-g|To|Pi-q)*. 



(68) 



The latter can be evaluated by means of Eq. (|66p . For 
q = we have immediately 



S.{Pf,Pi;0) = 6 



Pf-Pf 



\f{Pf^Pt)\' 

C^tot iP^) \P^ 



(69) 



while for q ^ Q an. extension of the rule (l66|) to differ- 
ent pairs of incoming and outgoing relative momenta is 
required. It can be constructed by formally taking the 
square root of ([66|l . Insertion into ([68|) brings about 
the square root of a product of two energy conserving 

2 _ 2 

(5-functions with arguments ^ .^ ' ± [pf — Pj) • q. Like 
with the delta- functions in Sect. lVBl this product implies 
that the parallel component q|| of the momentum separa- 
tion must be zero, thus restricting a q-integration to the 
plane perpendicular to the momentum change pj — Pi, 



i {Pf,Pz;q) 



i{Pf,Pi-qi_) if 911 =0 



otherwise. 



(70) 



The formal square root of the product of delta func- 
tions then reduces to a single proper Dirac function 

^ ((p/ Z^)' ^^'^ obtain, as the natural gener- 

alization of (p^ . 



C {Pf,Pt;qi_) =s 



p)-p1\ f {pf + q±,Pi + qi_) 



2 / VcTtot {pi + q±) \pi + q± 
f* {Pf - q±^Pi - qi_) 



\/crtot {Pi - q±) \Pi - gj 



(71) 



Inserting this expression into ([S7|) . together with rates 
determined by (|44p. one arrives directly at the complex 
rate density given by Eqs. (|1T|1. (|^ . 

We emphasize again that, compared to the microscopic 
phase space description in the preceding section, the line 
of reasoning is here quite different, and indeed more 
heuristic. The fact that the two lines of argument yield 
identical results indicates that their specific assumptions 
do reflect the underlying physics. 



VII. THE GAS INDUCED ENERGY SHIFT 

We now turn to the second part of the master equation, 
given by the superoperator TZ defined in Eq. ([5]). This 
term describes the coherent modification of the tracer 
particle dynamics due to the presence of the gas. As with 
the incoherent part £ given in a naive evaluation 
would take the expressions (fTT|) and for the rate 
and scattering operators at face value, and would thus 



yield an ill-defined normalization involving a ^ (0) /|r2| 
term. The correct normalization will be obtained in this 
section by implementing the appropriate restriction to 
the incoming states in the same way as in Sect. |Vl 

It should be noted that the effect of the energy shift 
described by TZ can usually be neglected when the inco- 
herent effects of the master equation play a role so that 
C dominates the dynamics. However, one can set up 
atom interferometer experiments where one beam inter- 
acts with a gas filled region such that only those atoms 
contribute to the detected signal which did not change 
their momentum by a collision. In this case, the effect 
can be measured as a gas-induced phase shift, and it is 
usually accounted for by attributing a refractive index to 
the gas (HSIMHlsil. 

Before startiiig the calculation, let us note 
that in Ref. 9* a slightly different term was 
given for the coherent modification, namely 
TZ'p = iTrgas([Re(T),ri/2[p,55pgjri/2]). it dif- 

fers from (O in the location of one of the f^/^ operators. 
In fact, the two superoperators TZ and TZ' yield the same 
gas induced energy shift when applied to the immobile 
system discussed in . For the present case of a tracer 
particle, TZ' has the disadvantage of introducing a weak 
dependence on P/P' which would need to be removed 
as an additional approximation. It therefore seems more 
natural to start from the form ([5]) right away, which is 
manifestly unitary. 

The calculational procedure can be carried out in com- 
plete analogy to the reasoning in Sect. |Vl Inserting the 
stationary state ([S]) into the expression ([5]) for TZp yields 
immediately the coherent modification of the evolution 
due to the presence of the gas. In momentum represen- 
tation it takes the form of Eq. (pS)) with the energy shifts 
given by 



En (P) 



{2TThY 



dPo A* (Po) 



(72) 



X (rel (po, P) I Vl" Re (To) Vy'\ rel (po, P)) ■ 



We can again switch to the center-of-mass frame by in- 
troducing the relative momentum 



Pn 



rel iPo,P) 



(73) 



as a function of Pq. This way the energy shifts take the 
form of an average over the gas momentum distribution 
function fi, 



En{P) = J dpo A* (Po) en (p„) 



(74) 



Like in the incoherent case, the function to be averaged 
can again be written as an expectation value with re- 
spect to a normalized momentum state of the relative 
motion — |lp„))((p„||- The function has the unit of 
an energy. 



e„(p„) = -/iTr(rJ/'Re(To)ry'pp„ 



(75) 



13 



When evaluating the expectation value, the restriction 
to incoming wave packets can again be implemented by 
replacing pp^ with its restricted version p'^ . It is given 
by Eq. ([55]) with replaced by p„, see Fig. O One thus 
obtains 



X exp 1 1 
X S{p, 



,x ■ w 

h 

w) Re 



r, 



da;||p„ 




|ApJ ^ 




\Pn 


W 



da; 



-LP„ 



-.1/2 




Pn - TT 



/ [Pn 



,Pn 



2 Y 

Like in Sect. |Vl we identify the cylinder base area with 
the scattering cross section, |5]p^| = a{Pn), and we 
note that the delta- function removes the component of 
w which is parallel to p„, so that the dependence on 
a; Hp vanishes in the integrand, 



en ip„ 



gas 



TO 

dx 



dw S 



Pn 



(27rft) 
da; 



-Lp„ 

— I exp 



w 1 



Pn 



Re 



w 



v/|Pn + fl|P„-f I 



^ (P„ + f ) g (p„ - f ) 



Carrying out the a;j|p^-integration, dx\\ 



|ApJ, 



one observes that the remaining a;^p^ -integration yields 
an approximate two-dimensional delta-function in w±p . 
Combined with the delta function in p„ • w/pn = w\\p 
this gives a three-dimensional S (w), which permits to do 
the lo-integration. One thus obtains 



e„ (p„ 



-2T:h 



2 "gas 



Re[/ {Pn:P„ 



(76) 



It shows that the energy shift is essentially determined 
by the real part of the forward scattering amplitude, a 
fact that is well-known in the field of neutron and atom 
optics. Its effect is often expressed by introducing an 
index of refraction ni, as discussed in Sect. IIXEI 



The calculations in Sects . fVl and I VII both indicate that 
TOin, the rate function in the center of mass frame, is 
given by Eq. (|42| . The complex rate Min, which de- 
termines the incoherent evolution in momentum repre- 
sentation according to p7|) . is then obtained by averag- 
ing TOin with the gas momentum distribution function 
p. Specifically, Eq. jHH) tells that Mn {P, P' ; Q) = 
/ dpf^p (pg) TOin (PfiPijO) with the relative momentary, 
p^, and q defined in However, the resulting ex- 

pression is not of the factorized form (|27|) needed below 
when stating the master equation in its operator repre- 
sentation. 

To arrive at (P7|) we first change the integration 
variable from Pq to p^. Moreover, the p distribution 
can be split symmetrically into a product of square 
roots, p{Pq) = ^^/^ (Po) M^^^ (Po)' since pg can be 
equally expressed as p^ + (^pj-- + P) m/M + qm/m^, or 
as + {pj + P') m/M - qm/m^, see ([231) -dSS])- Noting 
I det (9po/<9pj) I = m^/ml we have 



{P,P';Q) 



TO 
TO, 

TO 

M 

X TOin {Pf,Pi;q± 



Pi 



X ^1/2 , 



(Pf + P' 



TO 

M 



{Pf^p) 



-<1± 



TO* 



where we replaced q by qj^ in the arguments of p^^^, 
in accordance with Eq. (|48p . Having implemented the 
property (|4ip of the generalized function TOin, we can 
now insert its explicit form (|42p . which introduces the 
scattering amplitudes and an energy conserving delta- 
function, 

Mi„(P,P';Q) 

m \ ^ f , 1 /9 / m , , TO 

^'^'U+^iPf+pi-^^^^ 



9 9 

Pf - pf 



X — / (p/ + q^,Pj + q^) /* (p/ - q±,p^ - q^ 

TO: / \ J 



(77) 



VIII. OPERATOR REPRESENTATION OF THE 
QUANTUM LINEAR BOLTZMANN EQUATION 

So far, the derivation of the master equation was dis- 
cussed in the momentum representation. Let us now 
turn to the question how to obtain the quantum lin- 
ear Boltzmann equation in a representation-independent 
form. The result will then immediately prove the com- 
plete positivity and the translational covariance of the 
dynamical map defined by the master equation. 



From here it is a small step to arrive at the explicit ex- 
pression given in Eq. (|24p . In order to obtain a factorized 
expression we rather perform another change of variables, 

TO TO P±Q + P'±Q TO Q 

Pi P^ — P^ + T7 o 

TO* Al 2 TO* 2 

Due to its dependence on the transverse P and P' com- 
ponents this transformation has the remarkable effect of 
producing an integrand which is a product of P- and 
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P'-dependent factors: 
Min (P, P'; Q) 



TO* ./ V M " V MJ 2 



m \ Q 



x^^^/^{p+^P\, + {l-^]:^] 6(—p-Q 



TO \ Q \ g 



x/(rel(p,Px)-y,rel(p,Px) + | 
xr (rel(p,PY)-f,rcl(p,Pl) + | 



dp (5 



L{p,P Q;Q)L* {p,P'-Q;Q). 

(78) 



The second equality, which brings about the functions 
L {p, P;Q) defined in Eq. ([^5]) . emphasizes the factor- 
ization. Observing that the delta function restricts the 
p-integration to the plane = {peIR'^:p-(5 = 0} 
perpendicular to Q finally leads to the expression an- 
nounced in Eq. (|27p , since for any function g (p) 



dp S 



p Q 

Q 



dip) 



dpg{p). (79) 



The function L from Eq. ([28| is clearly well-suited to 
characterize the master equation, since it contains all the 
details of the coUisional interaction with the gas. It com- 
prises the elastic scattering amplitude / (p/,Pi) defined 
by the two-body interaction, the mass M of the tracer 
particle, the momentum distribution function /i (p) of the 
gas, its mass to and number density rigag. Unsurprisingly, 
the function L plays a central role for the operator rep- 
resentation of the master equation as well. It permits to 
define a family of jump operators acting in the Hilbert 
space of the tracer particle. 



-Q.p 



e^x-Q/fi^ (p, P; Q) , 



(80) 



where X and P are the corresponding position and 
momentum operators (and the function L is given in 
Eq. (|28pV The first factor effects a momentum exchange 
by Q, since exp (iX • Q/h) |P) = \P + Q), while the ap- 
pearance of P in the second factor renders the function 
L operator valued. This implies that both the scatter- 
ing amplitude and the momentum distribution function 
attain an operator character in (|28p . which is possible 
because the P-dependence of L will be analytic for any 
physically reasonable interaction potential. 

With the jump operators ((80| at hand it is straightfor- 
ward to construct the superoperator C, whose momen- 
tum representation is given by Eq. (|17p with Min from 

(123, 



£p = 



(81) 



This is the equation given in Ref. |10| (up to a trivial 
change of notation) . 

It is reassuring to observe that the form of the genera- 
tor (|8ip is in accordance with the most general structure 
of a translation-invariant and completely positive master 
equation as characterized by Holevo [32|, see [ssL l33 for 
a discussion. However, the summation in Ref. [32l | is 
here replaced by the p-integration over the plane Q"*" in 

m- 

A further consistency requirement is based on the 
transformation to a moving frame of reference. Denot- 
ing the velocity boost by V, the transformed state of the 
tracer particle is given by 



Pv 



ix-MV /h -ix-MV /n 
c pe , 



and the incoherent evolution in the new frame of refer- 
ence £v is thus related to C by 



^-iX-MV/h , ^iX-MV/h 



-iX-MV/h 



(82) 



However, the same super-operator must be obtained if 
we actively shift the momentum distribution fi (p) of 
the background gas, by setting pv (p) = p{p ~ mV) 
in the function L defining the jump operators ([50]) . The 
reason why this transformation of the gas motion must 
have the same effect as (|82|) is that the interaction be- 
tween the tracer particle and the gas depends only on 
their relative motion. Indeed, the functions L and Ly, 
based on the gas distributions /i and pv in ((28|) . are 
related by L{p,P- MV; Q) ^ Lv {p + niV±Q,P; Q). 
Noting also that a change of the integration variable 
p p'j^ ~ p + mV±Q in (|5T|) is possible, since it leaves 
the plane invariant, one easily proves the equivalence 
of the coordinate transformation and the shift of the mo- 
mentum distribution. 

As a final step, let us also incorporate the coher- 
ent modification of the tracer dynamics as discussed in 
Sect. IVIII The energy shift operator 



H, 



EniP) 



(83) 



is given by the operator- valued version of Eq. (jSOp . It 
permits to write the coherent modification part of the 
master equation ([29|) as Tip = {ifi)~^ [Hn, p] . This super- 
operator has the same invariance and transformation 
properties as discussed above in the case of C. In partic- 
ular, its transformation to a moving frame of reference 
analogous to ([82|) is equally obtained by replacing p with 
pv in (150)1 . 

To summarize this section we include the free mo- 
tion Hamiltonian H = P^/2M, thus writing the com- 
plete quantum linear Boltzmann equation ([3|) in the 
representation- independent form 



dtp =4 
in 



■ p2 

V- Hn,P 



2 / dQ 



dp{ 



-Qv 



(84) 



15 



IX. LIMITING FORMS 

As an important cross-check of the QLBE derived 
above, let us now see whether taking suitable limits re- 
duces its form to that of previously established equations. 



A. Classical linear Boltzmann equation 

The most obvious limiting motion is that of a classical 
particle. If all off-diagonal elements vanish in a motional 
state, {P\p\P') = 0, it is characterized by the diagonal 
momentum distribution /p (P) = {P\p\P) alone, and in- 
sofar it is indistinguishable from a classical state. One 
expects that the motion of the diagonal elements pre- 
dicted by ([5T|) is equal to the one described by the clas- 
sical linear Boltzmann equation. 

As follows from the discussion in Sect. IIIIl the QLBE 
implies that the diagonal elements /p (P) satisfy 



9r"/p(p) 



j dQMi{P;Q)MP-Q) 

^V(P)/p(P), (85) 



with the rates M^i^ (P) and Mj^] (P; Q) given by 
Eqs. (Uni) and ([22]). The notation df°^^ indicates that wc 
focus here only on the differential change in time which 
is due to the collision part C of the master equation. 

This equation should be compared to the classical lin- 
ear Boltzmann equation [2| for the momentum distribu- 
tion function fp{P). The traditional form of the collision 
integral reads, in our notation. 



dpdn 



rel (p, P) 



XCT (|rel (p, P)| n, rel (p, P)) 
x{m(p')/p' (P')-m(p)/p^'(P)} 



(86) 



where n is the unit vector of an angular integration with 
dn the associated element of solid angle. The values of 
P' and p' are determined by momentum conservation, 
granting in particular |rel(p',P') | = |rel(p, P) |. Us- 
ing the PT-invariance of the differential cross section, 
<^{PfiPi) = '^{PiiPf)i the classical linear Boltzmann 
equation can thus be rewritten in the explicit form 



dr''f;\p) 



n. 



'gas 



dpdrii I rel (p, P 

TO* ./ 

X 

xcr 



(87) 



^(p-rel (p, P) + |rel (p,P)|n,) 
a(rel(p,P),|rel(p,P)|n,) 
.:/p^' (P-frel(p,P)-|rel(p,P)|n,) 

-/p^'(P)^ f dpdn/ 1 rel (p,P) I 
x^(p)cr(|rel (p, P)|n/,rel (p,P)) . 



The angular integrations can be converted into three- 
dimensional integrals with a delta function. Noting 

|rel(p,P)|<5(|p,,/|-|rel(p,P)|) 



\p,j\'-\re\ip,P) 



one arrives, after the substitutions Pj j — > Pi-f — ' 
rel (p, P) — Pi J, at the form 

5rVp'(j') = y dP. Af ^1 (p,: ^ p) r; (p,) 

-/p^'(P) 1 dPsM^\P^Ps) 



with the classical rate density for the change of the tracer 
particle momentum from P,; to P f given by 

(P, ^ Pf) 

= — / dpo^l{Po)'J{Ye\{po,P^) + P^~Pf,re\{po,P^)) 



X 5 



|rel (po,P,;) + P, 



rel(Po,P.)r 



(89) 



It is now easy to see that the form (|88|) of the clas- 
sical linear Boltzmann equation is indeed identical to 
the diagonal part ^ of the QLBE, with M^\^ (P) and 
Mf„i (P; Q) given by Eqs. ([20 | and (|22l) . 



B. Pure collisional decoherence 

Another possible effect of the gas on a quantum tracer 
particle, and in a sense the other extreme compared to 
the classical dynamics on the diagonal, is the appearance 
of pure collisional decoherence. It follows from the QLBE 
([5T|) in the limit where the mass M of the tracer particle is 
much larger than the mass m of the gas molecules, so that 
there is no energy exchange during a collision. Taking 
m/M to zero simplifies the function (|28p characterizing 
the jump operators in (j80p . and renders it independent 
of P, 



L{p,P:Q) 



''gas 



Q 

2 



1/2 



xf[p±Q-^,p±Q + ^)-m 
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It follows that the generator of the incoherent evolution 
nil reduces to the form* 



TT^o n,. 



'9 9 

Pf-Pi 



<^{Pf:P^) 



This is the master equation of pure coUisional decoher- 
ence discussed by Gallis and Fleming and derived in 
its final form in Ref . @ . It describes an exponential decay 
of the off-diagonal elements in position representation, 

{X\Cp\X') -F{X-X'){X\p\X'), (92) 

with a localization rate given by 



F{R-R') = 



Pj-P) 



xa(p„p,){l-e<«-«')(p.-P/)A}. 

(93) 

This loss of coherence in the position basis can be at- 
tributed to the amount of position information (or 'which 
path' information) gained by the colliding gas. Recently, 
it has been observed that interfering fuUerene molecules 
display a reduction of interference visibility in agreement 
with this equation 0, . 



Specialization to the Maxwell-Boltzmann 
distribution 



So far, we kept the momentum distribution p of the 
gas molecules unspecified. This served to highlight the 
generality of the equations and it permitted, at the end 
of Sect. I Villi to discuss the implications of a transforma- 
tion of the momentum distribution. However, the most 
important choice is of course that of a Maxwell gas, char- 
acterized by a temperature T = l/f3k^. The remaining 
discussions of limiting forms in this section will be done 
with the corresponding Maxwell-Boltzmann distribution. 



M/3 (P) 




(94) 



where pp — \/2m/ [3 is the most probable momentum. 
We note that the statistical operator of the gas then takes 
the form 



(95) 



* In the same limit m/M — > the energy shift operator Hn from 
H83|l turns into a constant, so that it has no observable conse- 
quences for a constant gas density. 



with Ath = \/ '2tt1t? j3 / m the thermal de Broglie wave 
length and Iq the projectors to the normalization region, 
which are known from ([6|. 



D. Wealt coupling result 

A first limiting form of the QLBE that was obtained for 
the Maxwell-Boltzmann distribution is the weak coupling 
result by one of us [HI, [13, [I3l ■ Its derivation differs 
strongly from the approach of the present article, using 
the van Hove expression to relate the dynamic structure 
factor of the gas to the differential cross section in the 
laboratory frame. 

It can be regained from the present QLBE by replac- 
ing the exact scattering amplitude / in pS]) by its Born 
approximation /s, which is proportional to the Fourier 
transform of the interaction potential. 



fB{Pf,Pr) ^~4:Tr^hm4pj:\V (x) |p. 



dx V {x) exp 



(96) 



[Pf ^ P^ 



X 



Importantly, the Born approximation depends only on 
the momentum transfer — p^, but not on the energy 
in the relative motion. Even though /b violates the uni- 
tarity relation expressed by the optical theorem, it can 
be used to approximate the scattering amplitude if the 
energy of the relative motion is much larger than the 
interaction energy. 

Inserting the Born approximation (|96p into the 
function (|28p defining the jump operators Lq p — 
^iX Q/hj^ (p, P; Q) from Sect. lVml one notes that the P- 
dependence drops out in the scattering amplitude. As a 
result, the Born approximation of the function (j28p can 
be written as 



Lb{p,P;Q) = 



^gas 



1/2 



/b i-Q) VSiQ,P) 



exp 



pIq 



(97) 



where S (Q, P) is the dynamic structure factor of the 
Maxwell gas [33| , 



S{Q,P) 



I f3m 1 

exp 

27r Q ^ 



((1 



M 



9 m p 



QY 



8mQ^ 



(98) 



Since the p-dependence in ([971) appears just as a factor, 
one can carry out the Q -integration in the operator 
representation (j8ip of the master equation. The weak 
coupling approximation of the QLBE thus reduces to the 
form 



J dQ |lqpL^^ - ipLt^LQ - ll^QlQpj (99) 
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The corresponding jump operators 
Lq = e^xQ/'* 



n 
m 



-,1/2 



^S{Q,P)aB (Q) 



, (100) 



are determined by the cross section in Born approxima- 
tion, as (Q) = I/b (Q) P = I/b (-Q) rather than the 
individual scattering amphtudes. As a result, the mo- 
mentum operator P shows up with a particularly simple 
functional dependence, given by the dynamic structure 
factor Recently, the behavior of this equation was 

studied by means of a Monte Carlo simulation [s^ . 

The weak coupling form of the QLBE coincides with 
the expression derived earlier by one of us in Ref. [l2| 
(as can be seen if one combines the Eqs. (2) and (25) 
of that article, setting i {q) = —fs (q) /An^hm^, and z = 
nAjjj). This agreement is quite remarkable, given the very 
different type of argumentation in 12], and it serves to 
corroborate the validity of the present result. 

Incidentally, Eq. ((99|) also shows that the full QLBE 
cannot be obtained from the weak-coupling result by 
simply replacing fs by the proper scattering amplitude. 
This procedure would be ambiguous since the exact scat- 
tering amplitude is not just a function of the momentum 
transfer. As discussed at the end of Sect. |Vl the depen- 
dence of the scattering amplitudes on the tracer particle 
momentum, which dropped out in the Born approxima- 
tion, is required if one wants to cover the full set of pairs 
of scattering trajectories allowed under energy and mo- 
mentum conservation. 

In this context, it is worth mentioning that similar 
equations are obtained from a heuristic method of dealing 
with products of delta functions like the ones encountered 
in Sec. IIVBI see e.g. @, HI]. There, one of the energy 
delta functions is replaced by a finite Fourier integration 
over the 'elapsed time' as is done in derivations of Fermi's 
golden rule. Effectively, this amounts to a treatment in 
second order perturbation theory where it is permissible 
to identify the interaction Hamiltonian with the Born ap- 
proximation of the T-matrix. Although this brings scat- 
tering theory language into the game, one should not be 
tempted to conclude that the non-perturbative equation 
can be obtained by using the exact T-matrix. 



E. Index of refraction 

An application of the QLBE involving a rather spe- 
cial limit concerns the interference of matter waves in a 
Mach-Zehnder setup, where two interference paths are 
spatially separated by a macroscopic distance. One may 
ask how the interference fringe pattern changes if the 
particle is allowed to interact with a background gas in 
one of the interferometer arms. This setup was realized 
experimentally with Na and Li atoms [U [2^ [sij (and 
it is a common configuration in neutron interferometry 
where the "background gas" consists rather of thermal- 
ized condensed matter pol l4l|). 



In these situations the beam is strongly coUimated, 
while the likelihood of double collisions is small, so that 
after any collision that changes the momentum of the 
interfering tracer particle the latter will be blocked by 
the interferometer apertures. As a consequence, only 
forward-scattered amplitudes may contribute to the in- 
terference pattern, thus making the energy shift (|30p di- 
rectly observable as a change in the phase of the inter- 
ference pattern. At the same time an attenuation of the 
recorded signal is observed. 

The phase shift is usually accounted for by attributing 
a real index of refraction ni to the gas, which describes 
the modification of the de Broglie wavelength due to the 
energy shift. Exploiting the analogy between the force- 
free Schrodinger equation and the Helmholtz wave equa- 
tion [13, [43] the index of refraction for matter waves is 
determined by the ratio of the energy shift En (P) from 
([50)1 to the vacuum kinetic energy E'kin = P^/2M of the 
particle. 



= 1 - 



1 

X Re 



Ekin {P) 



(101) 



/ 



1 



2m, 



dPo A* (Po) 



M(Po,J')]" 



Here we took a rotationally invariant scattering ampli- 
tude, f{pf,Pi) = f{6\E^ci), with 6 = <{pf,Pi) and 

The index of refraction is typically close to unity, and 
therefore well approximated by the linearization 



ni = l + 27r' 



ngas A/ 



Re(/), 



(102) 



where K = P/h is the wavenumber of the interfering 
particle and Re(/) denotes the real part of the thermally 
averaged forward scattering amplitude. 



(/> = 1 dpo M (po) / (^0; ^ [rel (po, P)?^ 



(103) 



It is common in optics to account for the absorption in 
a medium by introducing an imaginary part to the index 
of refraction, which describes the exponential decay of 
the beam intensity. In the case of a background gas the 
tracer particles do not get absorbed, of course. However, 
for a strongly coUimated particle beam one expects an 
exponential decay of the beam intensity after a distance 
L, since collisions with the background gas decrease the 
probability of remaining in the beam, thus reducing the 
fraction of particles taking part in the coherent, wave- 
like behavior. The decay may be described by neglecting 
the gain term in Eq. (|85p . and integrating the remain- 
ing equation d^°^^fp{P) = -M^i^ (P) /p(P) up to time 
t = L/V, with V — P/M, yields the reduction factor 
exp (— Af°"* (P) L/V). By comparing this to the damped 
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intensity of a wave, exp {—2n2KL), one finds 

M°i"* (hK) 



«2 



2hK^/M 



(104) 



Inserting M°"* from (PT|) and using the optical theorem 
[13], that is, PitTtot (Pi) = 4^^Im [/ {Pi,Pi)], one gets an 
expression analogous to (|102p . 



(105) 



with Im(/) the imaginary part of (|103p . It follows that 
the combined effect of the energy shift and the reduc- 
tion of the amplitude of the coherent beam can be de- 
scribed by a complex index of refraction n — rii + in2 = 



2^ngasM/ {m,K^) (/). 



In case of a Maxwell-Boltzmann distribution 
average takes the form 



the 



°° dvv . , 
sinh 

V0 y 



2vV 



exp 



x/(0;— «^ 



(106) 



with V = hK/M the velocity of the interfering particle 
and v/3 = pfj/m. This expression of the thermally av- 
eraged forward scattering amplitude coincides with the 
one obtained by C. Champenois and collaborators in 
Ref. [H, |3| with a very different argumentation. It is 
used in the analysis of the recent experiment with Li 



28| 
30| 



atoms [31|. We note that the earlier experiments [2 
and the corresponding theoretical treatments [2£ 
were based on different expressions for (/) which we con- 
sider incorrect, see also the discussion in 44]. 



F. Diffusive limit 

A final important border case is the diffusive limit 
which is applicable if the tracer state is close to ther- 
mal and if its mass is much greater than the gas parti- 
cle mass, so that the motion is characterized by small 
momentum transfers. As discussed in [45j, an expan- 
sion of the jump operators (|80p to second order in the 
tracer position and momentum operators X and P is then 
permissible. In the special case of a constant scattering 
cross-section = al°^^^/AT: the QLBE then 

transforms into the generalized Caldeira-Leggett master 
equation 



Cp 



'h2 



D 



i=l 
3 



(107) 



It differs form the original equation Ti'] in the presence of 
the last term on the r.h.s of (|107p . which is necessary to 



ensure the complete positivity of the dynamics generated 
by £ (46j . We emphasize that, unlike in derivations using 
phenomenological choices for the model environment |47l . 
S, m, 113] , the friction and diffusion coefficients ?/ and 
Dpp are here uniquely specified by physically measurable 
properties of the gas. Specifically, the calculation in [iH] 
shows that the friction coefficient rj is determined by the 
temperature, the mass, and the density of the g well 
as by the scattering cross section. 



3^1/2 



M 



(108) 



The momentum diffusion constant Dpp is related to rj by 
the fluctuation-dissipation relation 



(109) 



Moreover, the coefficient of the "position-diffusion" term 
Dxx, is already determined by rj and Dpp, and it is given 
by the smallest value compatible with complete positiv- 
ityS 0, 



h^(3 

16M 



{^]Dpp. 

\am) 



(110) 



This shows that the diffusive limit turns the QLBE into 
the closest possible quantum analogue to the correspond- 
ing classical Kramers equation [5l| . 



X. CONCLUSIONS 

As seen in the previous section, all relevant limiting 
cases of the QLBE lead naturally to established master 
equations. In conjunction with the detailed derivations 
presented in Sects. lIVHVIIi this provides ample evidence 
that Eq. ([84)1 is the appropriate full quantum analogue 
of the classical linear Boltzmann equation. As such, it 
serves to describe non-perturbatively and in a unified 
framework the effects of decoherence and dissipation on 
a tracer particle. 

One reason that seems to have prevented this equation 
from being formulated earlier is the curious appearance 
of a second momentum integral in (|8ip which, in addi- 
tion to the integration over the momentum exchange Q, 
runs over the plane perpendicular to Q. This makes the 
equation a bit cumbersome at first sight, at least if rep- 
resented in a specific basis. However, we have seen in 
the course of the derivation that this five-dimensional in- 
tegration is necessary if one wants to cover all the pairs 
of scattering trajectories which are allowed by both the 



Diosi's equation Q leads to the same structure (I107I I. but the 
"position diffusion" constant Dxx is a complicated function of 
the cross section instead of being simply related to Dpp and f). 
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energy and momentum conservation and by the choice of 
Q. From a quantum mechanical point of view it is indeed 
natural to expect that the full set of possible scattering 
amplitudes contributes to the dynamics. The somewhat 
unwieldy explicit form is then the inevitable result of the 
transformation from the center of mass frame, where the 
scattering transformation takes place, to the laboratory 
frame needed for the tracer particle coordinate. 

Needless to say, the QLBE has a number of limitations. 
Like the classical linear Boltzmann equation, it cannot 
be applied in environments where the central Markov as- 
sumption is inappropriate, such as liquids. Moreover, it is 
not applicable at temperatures where the gas is quantum 
degenerate, and it is far from obvious how this possibility 
could be incorporated in the framework of the monitoring 



approach. Finally, let us reiterate that we presented here 
a physical derivation which, though stringent and lead- 
ing to a uniquely distinguished equation, may be hard 
to substantiate from a formal point of view. An alterna- 
tive, mathematically more rigorous derivation would be 
certainly desirable. 
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